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Abstract 

Let a n = [1 + n • log 2 3J. For the Collatz 3x + 1 function exists for each 
n € N a set of different residue classes (mod 2°") of starting numbers s 
with finite stopping time ct(s) = a n . Let z n be the number of these residue 
classes for each n > 0 as listed in the OEIS as A100982. It is conjectured 
that for each n > 4 the value of z n is given by the formula 


Zn 


5(n —2) 
3 

n — 2 


n ~ 1 / I 3(n-i)+<5 I 

V n ~ i 


where 5 £ Z assumes different values within the sum at intervals of 5 or 
6 terms. This allows to create an iterative algorithm which generates z n 
for each n > 12. This has been proved for each n < 10000. The number 
210000 has 4527 digits. 
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1 The Collatz 3x + 1 conjecture 


The Collatz 3x + 1 function is defined as a function T : N —> N on the set of 
positive integers by 


T{x) 



if x is even, 

if x is odd. 


Let T°(s) = s and T k (s) = T (T k 1 (s)) for k € N. Then the Collatz sequence 
for s € N is C(s) = ( T k (s ) | k = 0,1, 2,3,...). 

For example, the starting number s = 11 generates the Collatz sequence 

C{ 11) = (11,17,26,13,20,10,5,8,4,2,1,2,1,2,1,...). 

A Collatz sequence can only assume two possible forms. Either it falls into a 
cycle or it grows to infinity. The unproved conjecture to this problem is that 
each Collatz sequence enters the cycle (1,2). 


2 Stopping time 

2.1 The stopping time cr(s) 

Collatz’s conjecture is equivalent to the conjecture that for each s € N, s > 1, 
there exists k € N such that T k (s) < s. The least k € N such that T k (s) < s 
is called the stopping time of s, which we will denote by cr(s). It is not hard to 
verify that 

cr(s) = 1 if s 
<7(s) =2 if s 
cr(s) =4 if s 
<j(s) =5 if s 
cr(s) =7 if s 
<j(s) =8 if s 
cr(s) = 10 if s 
and so forth. 

Let cr„ = [1 + n • log 2 3J then generally applies for all n € N, n > 0, that 
cr(s) = cr n if s = X\, X 2 , £ 3 , • • •, x z (mod 2 CT "). 

Let z n be the number of residue classes ( mod 2 <T ") for each n > 0. Then we have 
z 0 = 1, zi = 1, z 2 = 1, z 3 = 2, z 4 = 3, z 5 = 7, z 6 = 12, z 7 = 30, ... 


= 0 (mod 2), 

= 1 (mod 4), 

= 3 (mod 16), 

= 11,23 (mod 32), 

= 7,15,59 (mod 128), 

= 39,79,95,123,175,199,219 (mod 256), 

= 287,347,367,423,507, 575,583,... (mod 1024), 
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Theorem 1: Let a n = [1 + n ■ log 2 3J then for each n £ N, n > 0, there exists 
a set of z n residue classes ( mod 2 cr ’ 1 ) of starting numbers s with finite stopping 
time er(s) = a n . For each n > 6 there is z„ > 2 • z n -\. (End of theorem) 

Everett[1] proves that almost all k £ N have finite stopping time, and Ter¬ 
ras [3] gives a probability distribution function for stopping times. 

The possible stopping times er(s) are listed in the OEIS as A020914. The associ¬ 
ated residue classes ( mod 2 CT ") are listed in the OEIS as A177789. The number 
of residue classes z n for n > 1 are listed in the OEIS as A100982. 

Appendix 5.4 shows a list of the first residue classes ( mod 2 <T ’ 1 ) up to a(s) = 15. 


2.2 A stopping time term formula for odd s 

Theorem 2: Let C a (s) = (T k (s ) | k = 0,..., a) with a > 1 be a finite sub¬ 
sequence of C(s), and let a n = [1 + n ■ log 2 3j. Then for each ueNan odd 
starting number s has the stopping time cr(s) = cr n , if the appropriate subse¬ 
quence C r 7 "~ 1 (s ) consists of n odd terms, and cq = k, if and only if T k (s ) in 
C" T " _1 (s) is odd. Then, there is 


T^(s) 


3 n 


S “h 


i= 1 


< T°(s). 


(1) 


(End of theorem) 


Example: For n = 4 there is a 4 = [1 + 4 • log 2 3J = 7. For s = 59 we get 
by equation (1) 


T 7 (59) 


3 4 3 3 2 0 + 3 2 2 1 + 3 1 2 3 + 3°2 4 

59 + - V - 


= 38 < 59. 


Explanation: The subsequence C 6 (59) = (59,89,134, 67,101,152, 76) consists of 
four (n = 4) odd terms 59,89,67,101. The powers of two cq yield as follows: 
T° = 59 is odd, so a± = 0. T 1 = 89 is odd, so «2 = 1. T 2 = 134 is even. T 3 = 67 
is odd, so «3 = 3. T 4 = 101 is odd, so 0:4 = 4. T 5 = 152 is even. T 6 = 76 is even. 

Note: There is a(s) = 7 not only for s = 59, but for every s = 59 ( mod 2 7 ). 


3 Diophantine equations 

3.1 Subsequences and their binary simplification 

Theorem 3: Let m = j and a n = [1 + n ■ log 2 3J for each n £ N, n > 4. 

According to the conditions of theorem 2, if an odd starting number s has the 
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stopping time a(s ) = a n , then the first m + n terms in C(s ) represents suffi¬ 
ciently the stopping time of s, because all further terms are even before the term 
T rJn (s) < s is reached. (End of theorem) 

To simplify the distribution of the even and odd terms in C(s) let ”0” rep¬ 
resents an even term and ” 1 ” represents an odd term. 

Example: For n = 4 there is m = = 1 and 04 

Then for the subsequences C an (s) there is 

C 7 { 7) = (7,11,17,26,13,20,10,5) simplified by 

C 7 (15) = (15,23,35,53,80,40,20,10) simplified by 
C 7 (59) = (59,89,134,67,101,152,76,38) simplified by 

And for the sufficiently subsequences C m+n ~ 1 (s ) there is 
C 4 (7) = (7,11,17,26,13) simplified by 

C 4 (15) = (15,23,35,53,80) simplified by 

C 4 (59) = (59, 89,134, 67,101) simplified by 

3.2 Binary tuples 

Let m = f° r eac ^ 11 6 N, n > 4. 

Let A{n) be a binary (to + n — 2)-tuple for each n eN,n> 4, defined by 

A(n) (ai, , Umj 5 0 m+n—2)7 

with ai, • • • , a m := 0 and a m+1 , ■ ■ ■ , a m+n - 2 := 1 . 

Let A j(n) be the set of A(n) and all permutations in lexicographic ordering of 
A{n) for each n S N, n > 4, whereby j is the number of all these tuples for each 
n, calculated by 

= (rn + n- 2 )! = 

to! • (n — 2)! \ n — 2 J' 

Let B(n) be a binary (to + ?r)-tuple for each n GN, n > 4, defined by 

B{n) = (6i, • • • , b m+n ) with bi := 1, b 2 := 1. 

Let Bj(?i) be the set of j tuples B(n) for each n G N, n > 4, where 63, • • • , b m+n 
is equal to exactly one tuple of A j (n). 

Example: For n = 7 there is to = 3 and j = 56. Therefore we get the 8 -tuple 
A( 7) = (0,0,0,1,1,1,1,1). There are further 55 permutations in lexicographic 
ordering of .4(7). Because of this A 56 (7) contains 56 different 8 -tuples, and 
B 56 (7) contains 56 different 10-tuples B( 7). Appendix 5.2 shows these example 
in detail. 


(1,1,1,0,1), 
( 1 , 1 , 1 , 1 , 0 ), 
( 1 , 1 , 0 , 1 , 1 ). 


= [1 + 4 • log 2 3J = 7. 

(1,1,1,0,1,0,0,1), 
( 1 , 1 , 1 , 1 , 0 , 0 , 0 , 0 ), 
(1,1,0,1,LO,0,0). 
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3.3 The Diophantine equations and their solutions 

According to theorem 2, the behaviour of a Collatz sequence is clearly related 
to the way in which the powers of 2 are distributed among the powers of 3 in 
the term s n ~ z 2 ai of equation (1). 


According to theorem 3, for each n G N, n > 4, only for the j binary tuples 
of Rj (n) the conditions of theorem 2 and equation (1) are complied. 


Theorem 4: Let a n = [1 + n ■ log 2 3J. By interpreting the binary tuples of 
Bj(n) as such a binary simplification for the even and odd terms in C m+n ~ 1 (s), 
there exists for each n G N, n > 4, according to theorem 2, for each binary tuple 
of Bj (n) a Diophantine equation 



i=1 


( 2 ) 


which has exactly one integer solution ( x, y) for 0 < x < 2 an , so that y = 
T an (x) < x in C(x). (End of theorem) 


Theorem 5: For each n > 4 the j solutions ( x,y ) represent j different residue 
classes of starting numbers s with stopping time 4 < cr(s) < cr n . More precise: 
For each n > 4 the j solutions represent z n residue classes of starting numbers s 
with the stopping time er(s) = a n , and j — z n residue classes of starting numbers 
s with stopping time 4 < a(s) < a n . For the z n solutions ( x,y ) with a(x) = a n , 
there is x the smallest number of the residue class [x] 2 ^n. (End of theorem) 


(■ x,y )■ 


(1, 

1, 

0 , 

0 , 1 , 

1, 

1) 

(1, 

1, 

o , 

1, 

0 , 

1, 

1) 

(1, 

1, 

0 , 

1, 

1, 

o , 

1) 

(1, 

1, 

0 , 

1, 

1, 

1, 

0) 

(1, 

1, 

1, 

o , 

o , 

1, 

1) 

(1, 

1, 

1, 

0 , 1 , 

o , 

1) 

(1, 

1, 

1, 

o , 

1, 

1, 

0) 

(1, 

1, 

1, 

1, 

o , 

o , 

1) 

(1, 

1, 

1, 

1, 

o , 

1, 

0) 

(1, 

1, 

1, 

1, 

1, 

o , 

0) 


The next ta- 
and the appropriate integer solution 

( 211 , 202 ) 

(107, 103) 

(123, 118) 

(219, 209) 

(183, 175) 

(199, 190) 

( 39, 38) 

( 79, 76) 

(175, 167) 

( 95, 91) 


Example: For n = 5 there are j = 10 integer solutions ( x,y ). 
ble shows the 10 binary tuples of B 10 (5) 


Of these 10 solutions, there are z 5 = 7 solutions for which x has the stopping 
time cr(x) = = 8, and 10 — 7 = 3 solutions for which x has stopping time 

4 < a(x) < 8. 


It is 211 G [3]is, 107 G [ll]s 2 and 183 G [23]32- So these 10 solutions 
represent 10 different residue classes of starting numbers s with stopping time 
4 < tr(s) < 8. According to chapter 2, there is 

if s = 3 ( mod 16), 
if s = 11, 23 ( mod 32), 
if se 39,79,95,123,175,199,219 ( mod 256). 


cr(s) = 4 
er(s) = 5 
a(s) = 8 
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Note that for cr(s) = 8 the 7 residue classes are equal to the values of x. 

The next table shows the distribution of the solutions (x, y) to the possible 
stopping times cr(x) for n = 4,..., 12. 



n 

4 

5 

6 

7 

8 

9 

10 

11 

12 

o(x) 











4 


0 

1 

1 

6 

28 

36 

165 

715 

1001 

5 


0 

2 

2 

10 

42 

56 

240 

990 

1430 

7 


3 

0 

0 

3 

15 

18 

84 

360 

495 

8 



7 

0 

7 

28 

35 

147 

588 

840 

10 




12 

0 

12 

12 

60 

252 

336 

12 





30 

0 

0 

30 

150 

180 

13 






85 

0 

85 

340 

425 

15 







173 

0 

173 

173 

16 








476 

476 

476 

18 









961 

0 

20 










2652 

Sum 


3 

10 

15 

56 

210 

330 

1287 

5005 

8008 


How to read: "Sum” is equal to the value of j and means the sum of the values 
of cr(x) for each n. For n = 5 there are 1+2+7=10 integer solutions (x, y). From 
these solutions, there are 7 with <r(x) = 8, 2 with a(x) = 5 and 1 with cr(x) = 4. 
No entry is equal to ” 0”. 


The PARI/GP programs 1 - 3 in Appendix 5.1 make the algorithm of this 
chapter more clearly. Program 2 generates the table of page 6, and program 3 
computes the values of the upper table. Appendix 5.3 shows Bj(n) with integer 
solutions (x, y) for n = 4,..., 7. 


3.4 Analysis of the solutions 

When analysing the entries of the last table, we can find that for each er(x) the 
number of solutions for all n > 4 are integer multiples of z n . 



n 

4 

5 

6 

7 

8 

9 

10 

11 

12 

a(x) 











4 


0- 1 

1 ■ 1 

1 • 1 

6 ■ 1 

28 • 1 

36- 1 

165- 1 

715- 1 

1001 - 1 

5 


0-2 

1 • 2 

1 • 2 

5-2 

21 • 2 

28 • 2 

120 • 2 

495-2 

715-2 

7 


1 ■ 3 

0-3 

0-3 

1 • 3 

5-3 

6-3 

28-3 

120-3 

165-3 

8 



1 ■ 7 

0-7 

1 • 7 

4-7 

5-7 

21 • 7 

84-7 

120-7 

10 




1 ■ 12 

0 • 12 

1 • 12 

1 ■ 12 

5 ■ 12 

21 ■ 12 

28 ■ 12 

12 





1 • 30 

0-30 

0-30 

1 • 30 

5-30 

6-30 

13 






1 • 85 

0-85 

1 • 85 

4-85 

5-85 

15 







1 • 173 

0- 175 

1 ■ 173 

1 ■ 173 

16 








1-476 

1-476 

1-476 

18 









1 ■ 961 

0-961 

20 










1 - 2652 

Sum 


3 

10 

15 

56 

210 

330 

1287 

5005 

8008 


7 
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It is not hard to recognize that the multiplying factors of z n like 1,4, 5, 6,21,... 
are numbers of Pascal’s triangle. It seems that these numbers are given for each 
n by the following binomial coefficients for k £ N. Note the different colours. 

^ = 0,1,1,5,6,28,36,165,220,1001,1365,6188,8568,38760,... 

^ = 1,1,4,5,21,28,120,165,715,1001,4368,6188,27132,... 




^ = 1,3,4,15,21,84,120,495,715,3003,4368,18564,27132,... 

Now we can conjecture that for each n > 4 the value of z n is equal to the 
difference of j and the sum of n — 2 terms, where each term is given by the 
product of a special binomial coefficient and a number z n for n = 2,...,n — 1. 
Let 



where S £ Z assumes different values within the sum at intervals of 5 or 6 terms. 
The example for n = 13 on page 9 and the algorithm on page 10 will make this 
clear. 



The next table shows the conjectured values for n = 13,..., 19. These values 
are not computed with program 3, but generated with program 4. 



n 

13 

14 

15 

16 

17 

18 

19 

a(x) 









4 


4368 ■ 1 

18564 ■ 1 

27132 ■ 1 

116280 ■ 1 

490314 ■ 1 

735471 ■ 1 

3124550 ■ 1 

5 


3003 • 2 

12376 • 2 

18564 • 2 

77520 • 2 

319770 • 2 

490314 ■ 2 

2042975 • 2 

7 


715-3 

3003 ■ 3 

4368 • 3 

18564 ■ 3 

77520 ■ 3 

116280-3 

490314 - 3 

8 


495-7 

2002 ■ 7 

3003 • 7 

12376 ■ 7 

50388 ■ 7 

77520 ■ 7 

319770 ■ 7 

10 


120- 12 

495 ■ 12 

715 • 12 

3003 ■ 12 

12376 • 12 

18564 • 12 

77520 ■ 12 

12 


28-30 

120 • 30 

165 ■ 30 

715 - 30 

3003 ■ 30 

4368 - 30 

18564 ■ 30 

13 


21 • 85 

84-85 

120 • 85 

495 • 85 

2002 • 85 

3003 • 85 

12376 ■ 85 

15 


5 ■ 173 

21• 173 

28 • 173 

120 • 173 

495 ■ 173 

715 ■ 173 

3003• 173 

16 


4-476 

15 • 476 

21 • 476 

84 ■ 476 

330 ■ 476 

495 • 476 

2002 • 476 

18 


1 • 961 

4-961 

5-961 

21 ■ 961 

84 • 961 

120 • 961 

495 • 961 

20 


0•2652 

1 • 2652 

1 • 2652 

5 • 2652 

21 • 2652 

28 • 2652 

120 • 2652 

21 


1 • 8045 

1 ■ 8045 

1 •8045 

4 • 8045 

15 • 8045 

21 • 8045 

84 • 8045 

23 



1 ■ 17637 

0■17637 

1 ■ 17637 

4 ■ 17637 

5■ 17637 

21 ■ 17637 

24 




1 ■ 51033 

1 • 51033 

3 • 51033 

4 • 51033 

15 • 51033 

26 





1 • 108950 

1 ■108950 

1 ■ 108950 

4 ■ 108950 

27 






1 ■ 312455 

1 • 312455 

3 ■ 312455 

29 







1 ■ 663535 

1 ■ 663535 

31 








1 - 1900470 

Sum 


31824 

125970 

203490 

817190 

3268760 

5311735 

21474180 
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For n = 14,..., 19 we need two new binomial coefficients. Note the new colours. 


L^J 


= 2,3,10,15,56,84,330,495,2002,3003,12376,18564,77520, 


m 


= 2,6,10,35,56,210,330,1287,2002,8008,12376,50388, 77520,... 


Example: This example shows the working of equation (3) for n = 13. 



zis = 31824 - 4368 • 1 - 3003 • 2 - 715 • 3 - 495 • 7 - 120 • 12 - 28 • 30 
- 21 • 85 - 5 • 173 - 4 • 476 - 1 • 961 - 0 • 2652 

zi3 = 31824 - 17424 - 6355 = 8045. 


3.5 An iterative stopping time algorithm 

The results of the analysis of the solutions and equation (3) enables us to de¬ 
vise an iterative algorithm which generates with only twelve initial numbers 
Zj,, z \2 each further number z n and also the number of all integer solutions 
(x. y ) for 0 < x < 2’ 7 ’ 1 from the Diophantine equation (2) with same stopping 
time u(x) as seen in the tables on page 7 and 8. 

The next PARI/GP program 4 shows this algorithm, which outputs a list of 
the numbers z n for 12 < n < limit, as listed in the OEIS as simple list of 
A100982. 
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The correctness of this algorithm for the values of z n has been proved with a 
counting-algorithm by Roosendaal/Noe[ 4] for each n < 10000. The num¬ 
ber 2:10000 has 4527 digits. Maybe the algorithm must slightly be modified for 
higher values of n. But the basic structure as seen in equation (3) will never 
been changed. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 


{/* Program 4 - The Winkler stopping time algorithm */ 
limit = 10000; 

zn=[] ; for(k=l, limit, zn=concat(zn,k)); 

/* 12 initial numbers */ 

zn[l]=l; zn[2]=l; zn[3]=2; zn [4] =3; 

zn [5] =7; zn[6]=12; zn[7]=30; zn[8]=85; 

zn[9]=173; zn[10]=476; zn[ll]=961; zn[12]=2652; 

/* main algorithm */ 
for (n = 12, limit , 

Sum = 0; a= -1; 

for(k=0, n-11, a=a +(-1)~floor (2*(k-2)/3)); 

for(i=2, 6, Sum = Sum + binomial(floor((3*(n-i)+ a)/2) ,n-i)*zn [i]); 
for(i=7, 11, Sum = Sum + binomial(floor ( (3*(n-i) + a-l)/2) ,n-i)*zn[i]) 
for(h=0, floor((limit-12)/6), 

v = 12; w = 5; 

if ((h<58) && (frac ( (h-6)/9)= = 0) , w = 4); 
for(k=0, 5, if (h>=k*9+7, v—)); 
c2 =58; 

for(k=l, floor(n/20), 

if (frac(k/2) !=0 , d = 52; e=5, d = 61; e = 6) ; 
c1 = c2; c2 = cl + d; 

if ((h>=cl) && (h<c2) && (frac((h-6 + k*2)/9)= = 0) , w=4) ; 
for(f=0, e, if (h>=f*9+cl+l, v--)); 

); 

p = 0; q = 0; 

for(i=v+h*6, v+h*6+w, 
delta=a-2-h; 

Sum=Sum+binomial(floor((3*(n-i)+delta)/2),n-i)*zn[i]; 
z_n = binomial(floor (5*(n-2)/3) ,n-2)-Sum; 
r1=0; r2=0; p + +; 

for(t=0, floor(n/15), 
c 1 =5 ; c2 =6 ; d = 0 ; 

if(t>cl, rl=2); if(t>c2, r2=2); 
for(k=2, floor(n/50), 

if (frac(k/2)==0, d=6, d = 7); 
cl=cl+d; c2=c2+d; 

if (t>cl , rl=k*2); if(t>c2, r2 = k*2); 

) ; 

for (k=t*9-r2, t*9+8-rl, 

if((n>k*6-t) && (h==k-2) && (p==n-k*6+t), 
print(n" "z_n); zn[n]=z_n; q=l; break; 

); 

); 

if (q==1, break ); 

) ; 

if (q==1, break ); 

); 

if (q==1, break ); 

) ; 

) ; 

> 
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The lines 12, 13 and 28, 29 of program 4 show the main calculations of equation 
(3). The only complicated thing is the computation of the pattern of the lower 
and upper bounds of the index of summation i and the appropriate value of 5. 

The upper limits of the for loops in the lines 19, 31 and 34 are dependent 
on the value of n or limit and must be large enough. The bigger the values of 
these upper limits, the longer the runtime of the algorithm. For example, for 
the small fixed upper limits 2, 24 and 2 the algorithm computes up to n = 1024 
much faster. 

The PARI/GP programs 5 and 6 in Appendix 5.1 show this algorithm with 
different outputs as seen in the tables on page 7 and 8. 

For another algorithm which generates the values of z n by a different way see 
Winkler[7] and Winkler[10]. But this older algorithm from 2011 is not yet 
immature. 


4 Conclusion and some answers 

4.1 A general stopping time theorem 

We see that the fact, if an odd starting number s has finite stopping time is 
only dependent on how the first m even and n odd terms are distributed in 
C : m+n- i( s ) anc j C CTti ~ 1 {s). Now we are able to formulate a general stopping 
time theorem without the use of equation (1). 

Theorem 6: Let m = |^ 2 ( ra ~ 2 ) j anc [ <j n = [1 + n • log 2 3j. For all n G N, n > 4, 
an odd starting number s has stopping time 4 < cr(s) < <j n , if the binary simpli¬ 
fication of the even and odd terms in the subsequence C m+n ~ 1 (s) is equal to a 
binary tuple of Bj(n) and if the subsequence C cr "-~ 1 (s ) consists of n odd terms. 
(End of theorem) 


4.2 Why the Collatz conjecture is a Diophantine equation problem? 

According to theorem 1, the Collatz conjecture is true, if the set of the residue 
classes (mod 2 <T ") of starting numbers for all n > 0 is equal to N. According 
to theorem 4 and 5, all possible residue classes for the stopping times are given 
by the Diophantine equations as their integer solutions. Therefore the Collatz 
conjecture is true, if there exists for each necessary residue class an appropriate 
Diophantine equation. 
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4.3 What are the consequences, if there exists another cycle than (1,2) 
or a sequence with infinite growth? 

Then, according to theorem 6, there exist an odd starting number s without stop¬ 
ping time. This starting number s has the property that for each n £ N, n > 4, 
there is no accordance of the binary simplification of the even and odd terms 
in the subsequence C m+n ~ 1 {s) and a binary tuple of Bj(n), or the subsequence 
C" T " _1 (s) consists not of n odd terms. 

4.4 Why has a small starting number like s = 27 a so comparatively large 
stopping time cr(27) = 59? 

Because out of the finite sets of Diophantine equations for each n > 4 not until 
n = 37 a Diophantine equation has the solution x = 27 or the residue class 
27 ( mod 2 59 ). Another answer: Not until n = 37 the binary simplification of the 
even and odd terms in the subsequence C m+n ~ 1 (27) is equal to a binary tuple 
of Bj(n) while the subsequence C <Tr * _1 (27) consists of n odd terms. 


4.5 Future work on the Collatz graph 

My future work is about the self-reference of the Collatz graph. I have devised 
two different directed graphs 1 , which illustrate that there exists a self-reference 
of odd subsequences C a (s ), clarified by the same background colour between 
both graphs and arrows in the first graph. 

This means on the one hand that odd subsequences C'°(s) with a small number 
of steps to reach 1 are connected with subsequences £7° (2s ± 1) or C a (4s ± 3) 
with a big number of steps to reach 1. For example C 3 (15) = (15,23,35,53) and 
C 3 (31) = (31,47,71,107). 

And on the other hand shows the first graph why and where to find numbers 
in the second graph. For example shows the green arrow in the first graph why 
C 1 (27) = (27,41) must hit the subsequence C 4 (31) = (31,47,71,107,161). In 
the second graph the subsequence C 6 (27) = (27,41,31,47,71,107,161) is lo¬ 
cated in column 9. 

This is only a very brief and insufficient exposition of the theory of self-reference. 
For details of the construction of the second graph and its correctness see WiN- 
KLEr[6]. Click here for a picture of the second graph which contains all odd 
numbers up to 341, and here for the complete structure up to column 6. 


1 See figure 1 on the next page. 
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Figure 1: Two directed Collatz graphs 
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5 Appendix 

5.1 Algorithms in PARI/GP 

The programs 1-3 use the function ’’NextPermutation(a)”, which generates all 
permutations in lexicographic ordering of a binary tuple A(n). 

1 NextPermutation(a)= 

2 { 

3 i=#a-l; 

4 while(!(i<l || a [i]<a[i + 1]) , i--); 

5 if(i<l, return(O)); 

6 k=#a; 

7 while (! (a[k]>a[i]) , k--) ; 

8 t=a[k]; a[k]=a[i]; a[i]=t; 

9 for (k=i+l, (#a+i)/2, 

10 t = a[k]; a [k]=a[#a + l + i-k] ; a [#a + l + i-k]=t; 

11 ); 

12 return (a); 

13 > 


Program 1 computes all integer solutions ( x , y) for 0 < x < from the Dio- 
phantine equation (2) for n = 7. 

1 {/* Program 1 */ 

2 n = 7; 

3 m=floor(2*(n-2)/3) ; 

4 Sigma=floor(l+n*log(3)/log(2)) ; 

5 j = (m+(n-2))!/(m!*(n-2) ! ) ; 

6 /* vectors for binary tuple */ 

7 Alpha=[0,0,0,0,0,0,0] ; 

8 A=[0,0,0,1 ,1 ,1 ,1 ,1] ; 

9 B= [1,1,0,0,0,0,0,0,0,0]; 

10 z = 0; q = 0; 

11 until(z==j, z++; 

12 /* generate B(n) from A(n) */ 

13 for(i = l, m+n-2, B [2+i]=A[i]); 

14 /* determine the n values for Alpha [i] */ 

15 i = l; for(k = l, m + n, if(B[k]= = l, Alpha [ i]=k-1; i+ + )); 

16 /* calculate integer solutions from Diophantine equation */ 

17 for (x=0, 2~Sigma, 

18 Sum=0; for(i=l, n, Sum=Sum+3~(n-i)*2"Alpha[i]); 

19 y= (3''n*x + Sum)/2" Sigma ; 

20 if (frac(y)==0, print(B" ("x", "y")"); q++); 

21 ) ; 

22 A = NextPermutat ion(A); 

23 ) ; 

24 print("There are "q" integer solutions (x,y) for n=7."); 

25 > 
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Program 2 computes all integer solutions (x, y) for 0 < x < 2 fT ’ 1 from the Dio- 
phantine equation (2) for 2 < n < limit. 

1 {/* Program 2 */ 

2 limit=100; 

3 for (n = 2 , limit , 

4 m=floor(2*(n-2)/3); 

5 Sigma=floor(l+n*log(3)/log(2)); 

6 j=(m+(n-2))!/(m!*(n-2)!); 

7 /* generate vectors for binary tuple */ 

8 Alpha=[]; for(i=l, n, Alpha=concat(Alpha,i)); 

9 for(i=l, n, Alpha[i]=0); 

10 A— [] ; for(i = l, m+n-2, A=concat(A , i )) ; 

11 for(i=l, m+n-2, if(i<=m, A[i]=0, A[i]=l)); 

12 B=[]; for(i=l, m+n, B=concat(B,i)); 

13 for(i=l, m+n, if(i<=2, B[i]=l, B[i]=0)); 

14 z =0; q = 0; 

15 until(z==j, z++; 

16 /* generate B(n) from A(n) */ 

17 for(i = l, m + n-2, B [2+i]=A [ i]); 

18 /* determine the n values for Alpha[i] */ 

19 i=1; for(k=l, m+n, if(B[k]==l, Alpha[i]=k-l; i++)); 

20 /* calculate integer solutions from Diophantine equation */ 

21 for(x=0, 2~Sigma, 

22 Sum=0; for(i=l, n, Sum=Sum+3~(n-i)*2"Alpha[i]); 

23 y= (S'*n*x + Sum )/2~ Sigma ; 

24 if (frac(y)= = 0, print(B" ("x", "y")"); q + + ); 

25 ) ; 

26 A = NextPermutat ion(A); 

27 ) ; 

28 print("There are "q" integer solutions (x,y) for n="n"."); 

29 print ; 

30 ) ; 

31 > 


Program 3 computes all integer solutions (x,y) for 0 < x < 2 an from the Dio¬ 
phantine equation (2) for 2 < n < 10 and counts the solutions with same stop¬ 
ping time g(x). 

1 {/* Program 3 */ 

2 for (n=2, 10, 

3 m=floor(2*(n-2)/3); 

4 Sigma=floor(l+n*log(3)/log(2)); 

5 j=(m+(n-2))!/(m!*(n-2)!); 

6 /* initialize stopping time counters */ 

7 wl = w2 = w3 = w4 = w5 = w6 = w7 = w8 = w9 = wl0=0; 

8 /* generate vectors for binary tuple */ 


9 

Alpha 

= []; fort 

i=l, n, Alpha=concat(Alpha, 

i)); 

10 


fort 

i=l, n, Alpha[i]=0); 


11 

A= [] ; 

for Ci = 1 , 

m+n-2, A=concat(A,i)); 


12 


for (i = 1 , 

m+n-2, if(i<=m, A[i]=0, A[ 

i]=1)) 

13 

B = [] ; 

for (i = l , 

m+n, B=concat(B,i)); 


14 


for (i = l , 

m+n, if(i<=2, B[i]=l, B[i] 

= 0)) ; 

15 

z =0; 

q = 0; 



16 

until(z==j, z++; 


17 

/ * 

generate 

B(n) from A(n) */ 


18 

f or (i = l, m+n 

-2, B [2+i]= A[i]); 


19 

/ * 

determine 

the n values for Alpha [i] 

*/ 
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20 

i=1; for (k=l, 

m + n , if (B[k]== 

= 1, Alpha [i]=k-l; i++)); 

21 

/* calculate integer solutions from Diophantine equation 

22 

for(x = 0, 2''Sigma, 



23 

Sum = 0; f or (i 

=1, n, Sum=Sum+3~(n- 

-i)*2"Alpha [i]); 

24 

y =(S^n+x+Sum)/2~Sigma; 



25 

if (frac(y)== 

0, q++; 



26 

/* calculate stopping time ST 

of x */ 

27 

ST = 0 ; c = x ; 




28 

unti1(c <x, 

if (frac(c/2)<>0, ST++; c = (3*c +1)/2) ; 

29 

if (frac(c/2)= = 0 , ST + + 

c=c/2; 

30 

if ((frac(c/2)<>0) && (c<x), break); 

31 

) ; 




32 

) ; 




33 

if ( ST = = 2 , 

wl + + ) ; if ( ST = = 

= 4, w2++); if (ST = = 5, w3 + + ); 

34 

if (ST = = 7, 

w4 + + ) ; if ( ST = = 

=8, w5++); if (ST==10, w6++); 

35 

if (ST = = 12, 

w7++); if(ST= 

==13, w8++); if (ST==15, w9++) 

36 

if (ST==16 , 

wlO + + ) ; 



37 

) ; 




38 

) ; 




39 

A = NextPermutat ion(A); 



40 

) ; 




41 

print; 




42 

print("There are 

"q" integer solutions (x,y) for n="n"."); 

43 

print("number of 

x with ST= 2 

"wl) 


44 

print("number of 

x with ST= 4 

" w2 ) 


45 

print("number of 

x with ST= 5 

" w3 ) 


46 

print("number of 

x with ST= 7 

" w4 ) 


47 

print("number of 

x with ST= 8 

" w5 ) 


48 

print("number of 

x with ST=10 

11 w6 ) 


49 

print("number of 

x with ST=12 

" w7 ) 


50 

print("number of 

x with ST=13 

" w8 ) 


51 

print("number of 

x with ST=15 

11 w9 ) 


52 

print("number of 

x with ST=16 

"wlO); 

53 

); 




54 

> 





Appendix 5.3 shows the output of program 2 for n = 4,..., 7. 
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Program 5 is equal to program 4 from chapter 3.5 but prints the number of all 
integer solutions as seen in the table on page 8. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 


{/* Program 5 - The Winkler stopping time algorithm */ 
limit=10000; zn=[]; for(k=l, limit, zn=concat(zn,k)); 

/* 12 initial numbers */ 

zn[l]=l; zn[2]=l; zn[3]=2; zn[4]=3; 

zn [5] =7; zn[6]=12; zn[7]=30; zn[8]=85; 

zn [9]=173; zn[10]=476; zn[ll]=961; zn[12]=2652; 

/* main algorithm */ 
for (n = 12, limit , 

Sum = 0; a=-1; 

for(k=0, n-11, a=a+(-1)~floor(2*(k-2)/3)); 

for(i=2, 6, Sum = Sum + binomial(floor ( (3*(n-i)+ a)/2) ,n-i)*zn [i] ; 
print(binomial(floor((3*(n-i)+a)/2),n-i)" * "zn [i]); 

) ; 

ford =7, 11, Sum=Sum+binomial(floor ((3*(n-i) + a-l)/2),n-i)*zn[i]; 
print(binomial(floor((3*(n-i)+a-l)/2) ,n-i)" * "zn [i]); 

) ; 

for(h=0, floor((limit-12)/6) , 
v = 12; w = 5; 

if ((h<58) && (frac((h-6)/9)==0), w=4); 
for(k=0, 5, if (h>=k*9+7, v—)); 
c2 =58; 

for(k=l, floor(n/20), 

if (frac(k/2) !=0 , d = 52; e=5, d = 61; e = 6) ; 
c1 = c2; c2 = cl + d; 

if ((h>=cl) && (h<c2) && (frac((h-6+k*2)/9)==0) , w = 4) ; 
for(f=0, e, if (h>=f*9+cl+l, v--)); 

); 

p = 0; q = 0; 

for(i=v+h*6, v+h*6+w, 
delta=a-2-h; 

Sum=Sum+binomial(floor((3*(n-i)+delta)/2),n-i)*zn[i]; 
if(n>12, print(binomial(floor((3*(n-i)+a-2-h)/2),n-i)" * "zn[i])) 
z_n=binomial(floor(5*(n-2)/3),n-2)-Sum; 
r1=0; r2=0; p + +; 

for(t=0, floor(n/15), 
c 1 =5 ; c2 =6 ; d = 0 ; 

if(t>cl, rl=2); if(t>c2, r2=2); 
for(k=2, floor(n/50), 

if (frac(k/2)==0, d=6, d = 7); 
cl=cl+d; c2=c2+d; 

if (t>cl , rl=k*2); if(t>c2, r2=k*2); 

) ; 

for (k=t*9-r2, t*9+8-rl, 

if((n>k*6-t) && (h==k-2) && (p==n-k*6+t), 
zn[n]=z_n; q=l; break; 

); 

); 

if (q==1, break ); 

); 

if (q==1, break ); 

); 

if (q==l, break ); 

) ; 

print("1 * "z_n); print; print; print; 

); 

> 
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Program 6 is equal to program 5 but prints the number of all integer solutions 
as seen in the upper table on page 7. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 


{/* Program 6 - The Winkler stopping time algorithm */ 
limit=10000; zn=[]; for(k=l, limit, zn=concat(zn,k)); 

/* 12 initial numbers */ 

zn[l]=l; zn[2]=l; zn[3]=2; zn[4]=3; 

zn [5] =7; zn[6]=12; zn[7]=30; zn[8]=85; 

zn[9]=173; zn[10]=476; zn[ll]=961; zn[12]=2652; 

/* main algorithm */ 
for (n = 12, limit , 

Sum = 0; a=-1; 

for(k=0, n-11, a=a+(-1)~floor(2*(k-2)/3)); 

for(i=2, 6, Sum = Sum + binomial(floor ( (3*(n-i)+ a)/2) ,n-i)*zn [i] ; 
print(binomial(floor((3*(n-i)+a)/2),n-i)*zn[i]); 

) ; 

ford =7, 11, Sum=Sum+binomial(floor ((3*(n-i) + a-l)/2),n-i)*zn[i]; 
print(binomial(floor((3*(n-i)+a-l)/2),n-i)*zn[i]); 

) ; 

for(h=0, floor((limit-12)/6) , 
v = 12; w = 5; 

if ((h<58) && (frac((h-6)/9)==0), w=4); 
for(k=0, 5, if (h>=k*9+7, v—)); 
c2 =58; 

for(k=l, floor(n/20), 

if (frac(k/2) !=0 , d = 52; e=5, d = 61; e = 6) ; 
c1 = c2; c2 = cl + d; 

if ((h>=cl) && (h<c2) && (frac((h-6+k*2)/9)==0) , w = 4) ; 
for(f=0, e, if (h>=f*9+c1+1, v--)); 

); 

p = 0; q = 0; 

for(i=v+h*6, v+h*6+w, 
delta=a-2-h; 

Sum=Sum+binomial(floor((3*(n-i)+delta)/2),n-i)*zn[i]; 
if (n>12, print(binomial(floor((3*(n-i)+a-2-h)/2) ,n-i)*zn [i])) 
z_n=binomial(floor(5*(n-2)/3),n-2)-Sum; 
r1=0; r2=0; p ++ ; 

for(t=0, floor(n/15), 
c 1 =5 ; c2 =6 ; d = 0 ; 

if(t>cl, rl=2); if(t>c2, r2=2); 
for(k=2, floor(n/50), 

if (frac(k/2)==0, d=6, d=7); 
cl=cl+d; c2=c2+d; 

if (t>cl , rl=k*2); if(t>c2, r2 = k*2); 

) ; 

for (k=t*9-r2, t*9+8-rl, 

if((n>k*6-t) && (h==k-2) && (p==n-k*6+t), 
zn[n]=z_n; q=l; break; 

); 

); 

if (q==1, break ); 

); 

if (q==1, break ); 

); 

if (q==l, break ); 

) ; 

print(z_n); print; print; print; 

); 

> 
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5.2 Permutations in lexicographic ordering for n — 7 

For n = 7 there is .4(7) = (0, 0,0,1,1,1,1,1), 

^ 56 ( 7 ) = 


o 

o 

o 

1 

1,1) , 

(1 

0 

1 

o 

o 

o 

o 

o 

1 

1,1) , 

(1 

0 

1 

o 

o 

o 

o 

o 

1 

1,1) , 

(1 

0 

1 

o 

o 

o 

o 

0 

1,1) , 

(1 

0 

1 

o 

o 

o 

o 

1 

0,1) , 

(1 

0 

1 

o 

o 

o 

o 

1 

1 ,0) , 

(1 

0 

1 

o 

o 

o 

o 

o 

1 

1,1) , 

(1 

0 

1 

o 

o 

o 

o 

o 

1 

1,1) , 

(1 

0 

1 

o 

o 

o 

o 

0 

1,1) , 

(1 

1 

0 

o 

o 

o 

o 

1 

0,1) , 

(1 

1 

0 

o 

o 

o 

o 

1 

1 ,0) , 

(1 

1 

0 

o 

o 

o 

o 

o 

1 

1,1) , 

(1 

1 

0 

o 

o 

o 

o 

0 

1,1) , 

(1 

1 

0 

o 

o 

o 

o 

1 

0,1) , 

(1 

1 

0 

o 

o 

o 

o 

1 

1 ,0) , 

(1 

1 

0 

o 

o 

o 

o 

0 

1,1) , 

(1 

1 

0 

o 

o 

o 

o 

1 

0,1) , 

(1 

1 

0 

o 

o 

o 

o 

1 

1 ,0) , 

(1 

1 

0 

o 

o 

(0,1 ,1 ,1 ,1 

0 

0,1) , 

(1 

1 

1 

o 

o 

o 

(0,1 ,1 ,1 ,1 

0 

1 ,0) , 

(1 

1 

1 

o 

o 

o 

(0,1 ,1 ,1 ,1 

1 

o 

o 

(1 

1 

1 

o 

o 

o 

o 

o 

o 

1 

1,1) , 

(1 

1 

1 

o 

o 

o 

o 

o 

o 

1 

1,1) , 

(1 

1 

1 

o 

o 

o 

o 

o 

0 

1,1) , 

(1 

1 

1 

o 

o 

o 

o 

o 

1 

0,1) , 

(1 

1 

1 

o 

o 

o 

o 

o 

1 

1 ,0) , 

(1 

1 

1 

o 

o 

o 

o 

o 

o 

1 

1,1) , 

(1 

1 

1 

o 

o 

o 

o 

o 

0 

1,1) , 

(1 

1 

1 

o 

o 

o 


®56 (7) = 


o 

o 

o 

1 

1,1) 


(1,1,1 

0 

1 ,0 , 1 , 1 

0,1) 

o 

o 

o 

1,1) 


(1,1,1 

0 

1 ,0 , 1 , 1 

1,0) 

o 

o 

o 

1 

1,1) 


(1,1,1 

0 

o 

o 

1,1) 

o 

o 

0 

1,1) 


(1,1,1 

0 

1 , 1,0 , 1 

0,1) 

o 

o 

1 

0,1) 


(1,1,1 

0 

1 , 1,0 , 1 

1,0) 

o 

o 

1 

1 ,0) 


(1,1,1 

0 

1 ,1 ,1,0 

0,1) 

o 

o 

o 

1 

1,1) 


(1,1,1 

0 

1 ,1 ,1,0 

1,0) 

o 

o 

o 

1 

1,1) 


(1,1,1 

0 

1 , 1 , 1 , 1 

o 

o 

o 

o 

0 

1,1) 


(1 , 1 , 1 , 1 

o 

o 

o 

1,1) 

o 

o 

o 


(1,1,1 

1 

o 

o 

o 

1,1) 

o 

o 

1 

1,0) 


(1,1,1 

1 

o 

o 

0,1) 

o 

o 

o 

1 

1,1) 


(1 , 1 , 1 , 1 

o 

o 

1,0) 

o 

o 

0 

1,1) 


(1,1,1 

1 

o 

o 

o 

1,1) 

o 

o 

1 

0,1) 


(1,1,1 

1 

o 

o 

0,1) 

o 

o 

1 

1,0) 


(1,1,1 

1 

o 

o 

1,0) 

o 

o 

0 

1,1) 


(1,1,1 

1 

o 

o 

0,1) 

o 

o 

1 

0,1) 


(1,1,1 

1 

o 

o 

1,0) 

o 

o 

1 

1,0) 


(1,1,1 

1 

0 , 1 , 1 , 1 

o 

o 

o 

o 

0,1) 


(1 , 1 , 1 , 1 

o 

o 

o 

1,1) 

(1,1,0,1,1,1,1 

0 

1,0) 


(1,1,1 

1 

o 

o 

0,1) 

(1,1,0,1,1,1,1 

1 

o 

o 


(1,1,1 

1 

o 

o 

1,0) 

o 

o 

o 

1 

1,1) 


(1 , 1 , 1 , 1 

o 

o 

0,1) 

o 

o 

o 

1,1,1) 


(1,1,1 

1 

o 

o 

1,0) 

o 

o 

0 

1,1) 


(1,1,1 

1 

1 , 0 , 1 , 1 

o 

o 

o 

o 

1 

0 , 1) 


(1,1,1 

1 

o 

o 

0,1) 

o 

o 

1 

1 ,0) 


(1,1,1 

1 

o 

o 

1,0) 

o 

o 

o 

1 

1,1) 


(1,1,1 

1 

1 , 1 , 0 , 1 

o 

o 

o 

o 

0 

1,1) 


(1,1,1 

1 

1 ,1 ,1,0 

o 

o 
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5.3 Bj(n) with integer solutions (x,y) for n — 4,... ,7 

This is also the output of program 2 for n = 4,..., 7. 


(1. 

1 , 

0 , 

1 , 

1) 


(59, 38) 




(1, 

1 , 

1 , 

o , 

1) 


(7, 

6) 





(1. 

1 , 

1 , 

1, 

0) 


(15 

, 10) 




There 

are 

3 

integer solutions (X 

, Y) for 

n = 4. 

(i. 

1 , 

0 , 

0 , 

1 , 

1 , 

i) 


(211, 202) 


(i, 

1 , 

0 , 

1 , 

0 , 

1 , 

i) 


(107, 103) 


a. 

1 , 

0 , 

1 , 

1 , 

o , 

i) 


(123, 118) 


(i, 

1 , 

0 , 

1 , 

1 , 

1, 

0) 


(219, 209) 


(i, 

1 , 

1 , 

o , 

0 , 

1, 

1) 


(183, 175) 


(i, 

1 , 

1 , 

o , 

1 , 

o , 

1) 


(199, 190) 


(i, 

1 , 

1 , 

o , 

1 , 

1, 

0) 


(39 

, 38) 



(i, 

1 , 

1 , 

1, 

0 , 

o , 

1) 


(79 

, 76) 



(i, 

1 , 

1 , 

1, 

0 , 

1, 

0) 


(175, 167) 


(i, 

1 , 

1 , 

1, 

1 , 

o , 

0) 


(95 

, 91) 



There 

are 

10 

integer 

solutions (X,Y) for n=5. 

(i. 

1 , 

0 , 

o , 

1 , 

1, 

1, 

i) 


(595 , 

425) 


(i, 

1 , 

0 , 

1, 

0 , 

1, 

1, 

i) 


(747 , 

533) 


(i, 

1 , 

0 , 

1, 

1 , 

o , 

1, 

i) 


(507 , 

362) 


(i, 

1 , 

0 , 

1, 

1 , 

1, 

o , 

i) 


(347 , 

248) 


(i. 

1 , 

0 , 

1, 

1 , 

1, 

1, 

0) 


(923 , 

658) 


(i. 

1 , 

1 , 

o , 

0 , 

1, 

1, 

1) 


(823 , 

587) 


(i, 

1 , 

1 , 

o , 

1 , 

o , 

1, 

1) 


(583 , 

416) 


(i. 

1 , 

1 , 

o , 

1 , 

1, 

o , 

1) 


(423 , 

302) 


a, 

1 , 

1 , 

o , 

1 , 

1, 

1, 

0) 


(999 , 

712) 


(i, 

1 , 

1 , 

1, 

0 , 

o , 

1, 

1) 


(975 , 

695) 


(i. 

1 , 

1 , 

1, 

0 , 

1, 

o , 

1) 


(815 , 

581) 


(i, 

1 , 

1 , 

1, 

0 , 

1, 

1, 

0) 


(367 , 

262) 


(i, 

1 , 

1 , 

1, 

1 , 

o , 

o , 

1) 


(735 , 

524) 


(i, 

1 , 

1 , 

1, 

1 , 

o , 

1, 

0) 


(287 , 

205) 


(i, 

1 , 

1 , 

1, 

1 , 

1, 

o , 

0) 


(575 , 

410) 


There 

are 

15 

integer 

solutions (X,Y) for n=6. 

a. 

1 , 

0 , 

0 , 

0 , 

1, 

1, 

1 , 

1 , 

i) 

(3523 , 

1883) 

(i, 

1 , 

0 , 

o , 

1 , 

o , 

1, 

1 , 

1 , 

i) 

(3827 , 

2045) 

(i. 

1 , 

0 , 

o , 

1 , 

1, 

o , 

1 , 

1 , 

i) 

(1299 , 

695) 

(i, 

1 , 

0 , 

o , 

1 , 

1, 

1, 

0 , 

1 , 

i) 

(979 , 

524) 

(i, 

1 , 

0 , 

o , 

1 , 

1, 

1, 

1 , 

0 , 

i) 

(2131 , 

1139) 

(i, 

1 , 

0 , 

o , 

1 , 

1, 

1, 

1 , 

1 , 

0) 

(2899 , 

1549) 

(i. 

1 , 

0 , 

1, 

0 , 

o , 

1, 

1 , 

1 , 

1) 

(3979 , 

2126) 

(i, 

1 , 

0 , 

1, 

0 , 

1, 

o , 

1 , 

1 , 

1) 

(1451 , 

776) 

(i, 

1 , 

0 , 

1, 

0 , 

1, 

1, 

0 , 

1 , 

1) 

(1131 , 

605) 

(i, 

1 , 

0 , 

1, 

0 , 

1, 

1, 

1 , 

0 , 

1) 

(2283 , 

1220) 

(i, 

1 , 

0 , 

1, 

0 , 

1, 

1, 

1 , 

1 , 

0) 

(3051 , 

1630) 

(i, 

1 , 

0 , 

1, 

1 , 

o , 

o , 

1 , 

1 , 

1) 

(187 , 

101) 

(i. 

1 , 

0 , 

1, 

1 , 

o , 

1, 

0 , 

1 , 

1) 

(3963 , 

2117) 

(i, 

1 , 

0 , 

1, 

1 , 

o , 

1, 

1 , 

0 , 

1) 

(1019 , 

545) 

(i, 

1 , 

0 , 

1, 

1 , 

o , 

1, 

1 , 

1 , 

0) 

(1787 , 

955) 

a, 

1 , 

0 , 

1, 

1 , 

1, 

o , 

0 , 

1 , 

1) 

(1755 , 

938) 

(i. 

1 , 

0 , 

1, 

1 , 

1, 

o , 

1 , 

0 , 

1) 

(2907 , 

1553) 

(i. 

1 , 

0 , 

1, 

1 , 

1, 

o , 

1 , 

1 , 

0) 

(3675 , 

1963) 

(i, 

1 , 

0 , 

1, 

1 , 

1, 

1, 

0 , 

0 , 

1) 

(1435 , 

767) 

(i, 

1 , 

0 , 

1, 

1 , 

1, 

1, 

0 , 

1 , 

0) 

(2203 , 

1177) 

(i, 

1 , 

0 , 

1, 

1 , 

1, 

1, 

1 , 

0 , 

0) 

(2587 , 

1382) 

(i. 

1 , 

1 , 

o , 

0 , 

o , 

1, 

1 , 

1 , 

1) 

(2007 , 

1073) 

(i. 

1 , 

1 , 

o , 

0 , 

1, 

o , 

1 , 

1 , 

1) 

(3575 , 

1910) 

(i, 

1 , 

1 , 

o , 

0 , 

1, 

1, 

0 , 

1 , 

1) 

(3255 , 

1739) 

(i, 

1 , 

1 , 

o , 

0 , 

1, 

1, 

1 , 

0 , 

1) 

(311 , 

167) 

(i. 

1 , 

1 , 

o , 

0 , 

1, 

1, 

1 , 

1 , 

0) 

(1079 , 

577) 

(i, 

1 , 

1 , 

o , 

1 , 

o , 

o , 

1 , 

1 , 

1) 

(2311 , 

1235) 

(i. 

1 , 

1 , 

o , 

1 , 

o , 

1, 

0 , 

1 , 

1) 

(1991 , 

1064) 

a, 

1 , 

1 , 

o , 

1 , 

o , 

1, 

1 , 

0 , 

1) 

(3143 , 

1679) 

(i. 

1 , 

1 , 

o , 

1 , 

o , 

1, 

1 , 

1 , 

0) 

(3911 , 

2089) 
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(1. 

1 , 

1 , 

0 , 

1 , 

1 , 

o , 

o , 

1 , 

1) 

(3879 , 

, 2072) 

(1. 

1 , 

1 , 

0 , 

1 , 

1 , 

o , 

1, 

0 , 

1) 

(935 , 

500) 

(1, 

1 , 

1 , 

0 , 

1 , 

1 , 

o , 

1, 

1 , 

0) 

(1703 , 

, 910) 

(1, 

1 , 

1 , 

0 , 

1 , 

1 , 

1, 

o , 

0 , 

1) 

(3559 , 

, 1901) 

(1. 

1 , 

1 , 

0 , 

1 , 

1 , 

1, 

o , 

1 , 

0) 

(231 , 

124) 

(1, 

1 , 

1 , 

0 , 

1 , 

1 , 

1, 

1, 

0 , 

0) 

(615 , 

329) 

(1. 

1 , 

1 , 

1 , 

0 , 

0 , 

o , 

1, 

1 , 

1) 

(3727 , 

1991) 

(1, 

1 , 

1 , 

1 , 

0 , 

0 , 

1, 

o , 

1 , 

1) 

(3407 , 

1820) 

(1. 

1 , 

1 , 

1 , 

0 , 

o , 

1, 

1, 

0 , 

1) 

(463 , 

248) 

(1, 

1 , 

1 , 

1 , 

0 , 

o , 

1, 

1, 

1 , 

0) 

(1231 

658) 

(1, 

1 , 

1 , 

1 , 

0 , 

1, 

o , 

o , 

1 , 

1) 

(1199 

641) 

(1, 

1 , 

1 , 

1 , 

0 , 

1, 

o , 

1, 

0 , 

1) 

(2351 

1256) 

(1, 

1 , 

1 , 

1 , 

0 , 

1, 

o , 

1, 

1 , 

0) 

(3119 

1666) 

(1. 

1 , 

1 , 

1 , 

0 , 

1, 

1, 

o , 

o , 

1) 

(879 , 

470) 

(1, 

1 , 

1 , 

1 , 

0 , 

1, 

1, 

o , 

1, 

0) 

(1647 

880) 

(1, 

1 , 

1 , 

1 , 

0 , 

1, 

1, 

1, 

o , 

0) 

(2031 

1085) 

(1. 

1 , 

1 , 

1 , 

1 , 

o , 

o , 

o , 

1, 

1) 

(2143 

1145) 

(1, 

1 , 

1 , 

1 , 

1 , 

o , 

o , 

1, 

o , 

1) 

(3295 

1760) 

(1, 

1 , 

1 , 

1 , 

1 , 

o , 

o , 

1, 

1, 

0) 

(4063 

2170) 

(1, 

1 , 

1 , 

1 , 

1 , 

o , 

1, 

o , 

o , 

1) 

(1823 

974) 

(1, 

1 , 

1 , 

1 , 

1 , 

o , 

1, 

o , 

1, 

0) 

(2591 

1384) 

(1, 

1 , 

1 , 

1 , 

1 , 

o , 

1, 

1, 

o , 

0) 

(2975 

1589) 

(1, 

1 , 

1 , 

1 , 

1 , 

1, 

o , 

o , 

o , 

1) 

(1087 

581) 

(1, 

1 , 

1 , 

1 , 

1 , 

1, 

o , 

o , 

1, 

0) 

(1855 

991) 

(1. 

1 , 

1 , 

1 , 

1 , 

1, 

o , 

1, 

o , 

0) 

(2239 

1196) 

(1, 

1 , 

1 , 

1 , 

1 , 

1, 

1, 

o , 

o , 

0) 

(383 , 

205) 

There 

are 

56 

integer 

solutions 

(X,Y) for n=7. 


5.4 Stopping time residue classes up to cr(s) = 15 

<r(s) = 1 

if s = 0 (mod 2) 

<t( S ) = 2 

if s = 1 (mod 4) 

<y(s) = 4 

if s = 3 (mod 16) 
o-(s) = 5 

if s E 11, 23 (mod 32) 

<t(s) = 7 

if sE 7, 15, 59 (mod 128) 
o-(s) = 8 

if sE 39, 79, 95, 123, 175, 199, 219 (mod 256) 
ct(s) = 10 

if sE 287, 347, 367, 423, 507, 575, 583, 735, 815, 923, 975, 999 (mod 1024) 
ct(s) = 12 

if s = 231, 383, 463, 615, 879, 935, 1019, 1087, 1231, 1435, 1647, 1703, 1787, 1823, 1855, 2031, 2203, 2239, 2351, 
2587, 2591, 2907, 2975, 3119, 3143, 3295, 3559, 3675, 3911, 4063 (mod 4096) 

cr(s) = 13 

if s E 191, 207, 255, 303, 539, 543, 623, 679, 719, 799, 1071, 1135, 1191, 1215, 1247, 1327, 1563, 1567, 1727, 1983, 

2015, 2075, 2079, 2095, 2271, 2331, 2431, 2607, 2663, 3039, 3067, 3135, 3455, 3483, 3551, 3687, 3835, 3903, 3967, 

4079, 4091, 4159, 4199, 4223, 4251, 4455, 4507, 4859, 4927, 4955, 5023, 5103, 5191, 5275, 5371, 5439, 5607, 5615, 

5723, 5787, 5871, 5959, 5979, 6047, 6215, 6375, 6559, 6607, 6631, 6747, 6815, 6983, 7023, 7079, 7259, 7375, 7399, 

7495, 7631, 7791, 7847, 7911, 7967, 8047, 8103 (mod 8192) 

cr(s) = 15 

if se 127, 411, 415, 831, 839, 1095, 1151, 1275, 1775, 1903, 2119, 2279, 2299, 2303, 2719, 2727, 2767, 2799, 2847, 
2983, 3163, 3303, 3611, 3743, 4007, 4031, 4187, 4287, 4655, 5231, 5311, 5599, 5631, 6175, 6255, 6503, 6759, 6783, 
6907, 7163, 7199, 7487, 7783, 8063, 8187, 8347, 8431, 8795, 9051, 9087, 9371, 9375, 9679, 9711, 9959, 10055, 10075, 
10655, 10735, 10863, 11079, 11119, 11567, 11679, 11807, 11943, 11967, 12063, 12143, 12511, 12543, 12571, 12827, 

12967, 13007, 13087, 13567, 13695, 13851, 14031, 14271, 14399, 14439, 14895, 15295, 15343, 15839, 15919, 16027, 

16123, 16287, 16743, 16863, 16871, 17147, 17727, 17735, 17767, 18011, 18639, 18751, 18895, 19035, 19199, 19623, 

19919, 20079, 20199, 20507, 20527, 20783, 20927, 21023, 21103, 21223, 21471, 21727, 21807, 22047, 22207, 22655, 

22751, 22811, 22911, 22939, 23231, 23359, 23399, 23615, 23803, 23835, 23935, 24303, 24559, 24639, 24647, 24679, 

25247, 25503, 25583, 25691, 25703, 25831, 26087, 26267, 26527, 26535, 27111, 27291, 27759, 27839, 27855, 27975, 

28703, 28879, 28999, 29467, 29743, 29863, 30311, 30591, 30687, 30715, 30747, 30767, 30887, 31711, 31771, 31899, 

32155, 32239, 32575, 32603 (mod 32768) 

and so forth. 
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